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M. Eschrig a , J. Ferrer h , and M. Fogelstr6m c 
a Materials Science Division, Argonne National Laboratory, Argonne, Illinois 60439 
b Departamento de Fisica, Facultad de Ciencias, Universidad de Oviedo, E-33007 Oviedo, Spain 
c Institute for Theoretical Physics, Chalmers University of Technology and Goteborg University, S-41296, Goteborg, Sweden 

(February 1, 2008) 

We show that in Sr2Ru04 the Fermi surface geometry as inferred from angle resolved photoemis- 
sion experiments has important implications for a pairing interaction dominated by incommensurate, 
strongly anisotropic, spin fluctuations. For a spin fluctuation spectrum consistent with inelastic neu- 
tron scattering experiments the system is close to an accidental degeneracy between even parity spin 
singlet and odd parity spin triplet channels. This opens the possibility of a mixed parity order pa- 
rameter state in S^RuO^ We determine the stable and metastable order parameter phases at low 
temperatures and discuss especially phases with order parameter nodes. 



The study of the superconducting material Si'2Ru04 
jl| has attracted considerable experimental and theoreti- 
cal attention because of its peculiar low energy electronic 
properties g] . It is a rare example of a multiband super- 
conductor with possibly spin triplet pairing symmetry 
H . Based on muon spin relaxation measurements jj] an 
early suggestion for the order parameter symmetry was 
the time reversal symmetry breaking s p x ±ip v '-state 
Recent experiments [pH7| on the contrary favor order pa- 
rameters with nodes pj . The spin excitation spectrum in 
Sr2RuC>4 as studied by inelastic neutron scattering (INS) 
reveals spin fluctuations in a crossover regime between 
ferromagnetic and anti-ferromagnetic fluctuations. The 
dominant contributions to the spin susceptibility, x(l)i 
are located around Q — (1 ± 6, 1 ± #, 0) 7r/a, with an 
incommensurability S ~ 0.4 |J. Similar results have 
also been suggested theoretically finding 5 « 1/3 [jl0| . It 
was suggested that in the cross-over regime between fer- 
romagnetic and anti-ferromagnetic spin fluctuations the 
order parameter changes from p-wave to d-wave |Io| , pd| . 
As experimental results become more reliable, it becomes 
possible to test if a pairing mechanism via incommensu- 
rate spin-fluctuations is consistent with a) the spin fluc- 
tuation spectrum as measured in INS, b) the electronic 
structure as tested by angle resolved photoemission and 
de Haas- van Alphen spectroscopy, and c) the experimen- 
tal implications for the order parameter symmetry. 

In this Letter we study the possible pairing symmetries 
for the order parameter in Sr2Ru04 resulting from a pair- 
ing interaction dominated by incommensurate, strongly 
anisotropic, spin fluctuations. Given the measured Fermi 
surface geometry in Sr2Ru04 JlJ,|l3|], we find that the 
system is close to two accidental degeneracies between 
the triplet E^ u superconducting state and one of the sin- 
glet superconducting states either in the B\ g (predom- 
inantly d-wave) channel or in the Ai g (predominantly 
g-wave) channel. This opens the possibility of a super- 
conducting state with mixed parity order parameter in 
Sr2RuC>4. The possibility of such superconducting states 
was discussed some time ago in connection with UPt3 
Q. We also study the different possible stable and 



metastable low-temperature phases in Sr2Ru04 for the 
three cases that either of the three relevant order param- 
eter symmetries E2 U , B\ g , or A\ g is dominant. We find 
both nodeless solutions and order parameters with nodes. 
A dominant triplet order parameter is supported only in 
a surprisingly small incommensurability region near the 
experimentally determined value S ~ 0.4. 

NMR experiments suggest a strong anisotropic spin 
susceptibility, with only the \zz peaked around the in- 
commensurate wave vector Q jf5[. Thus, we consider an 
anisotropic model with \ z = Xzz > Xxx = Xyy = X ± 
and neglect off-diagonal components (\xy = 0). For \ z 
we assume the following form (a = 1 in our units) 



X z (q)= £ 



XQ/4 
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2 SmI 



2x\ 
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This model susceptibility has three parameters. The 
overall magnitude, xq, determines the coupling strength 
in the dominant pairing channel, and thus the super- 
conducting transition temperature T c q. The other two 
parameters, the spin-spin correlation length, £ s a, and the 
degree of incommensuration, 5, determine the relative 
coupling strengths in the different symmetry channels, 
defining the symmetry of dominant and sub-dominant 
components. The specific form of the susceptibility al- 
lows a smooth cross-over from anti-ferromagnetic fluctu- 
ations, Qap = (1,1,0) 7r/a, to ferromagnetic, Qpu = 
(0,0,0), by tuning <5 from to 1. Extracting the values 
of £ s fl and 5 from the INS data || gives £ s a ~ 4.0a and 
5 w 0.4. The effective pairing interaction via spin fluc- 
tuation exchange is determined by the coupling function 
9(p)x l (P ~ p')9(p')- The coupling between spin fluctu- 
ations and quasiparticles, g{p), is approximated in what 
follows by a constant, g. To reproduce the structure 
of the experimentally probed three-sheet Fermi surface 
fl2| , [l3f we use the tight-binding dispersions 

€p = 2t l x COS Pa + 2ty cospy — 4i'' 1 cosp x cospy — p} . (2) 

The band index i labels the xy, xz and yz bands. The pa- 
rameters of the dispersions (t x , t y , t'' 1 , //) are taken from 
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Refs. [ plj| . Additionally, a hybridization of the xz and 
the yz bands is given by t± —0.1 eV Q . 

The order parameter in weak coupling BCS theory 
is determined by the well-known non-linear BCS self- 
consistency equation. The smallness of the gap in 
Sr2Ru04 of about 1 mcV [p| allows us to restrict mo- 
mentum summations in this gap equation to the vicin- 
ity of the Fermi surface. In a standard way this pro- 
cedure leads to a replacement of the momentum sum 
by a Fermi surface average, (• • -)p f , weighted with the 
angle resolved density of states, N(pj) = Nfi%(pf) = 
Nf \v f (p f )\- 1 /{\vf(p' f )\- 1 )p'. N f is the total density of 
states at the Fermi level and Vf is the Fermi velocity. We 
obtain Nf = 1.78 states/eV per spin and unit cell. Intro- 
ducing the dimensionless coupling function x — Nf g 2 x> 
the weak coupling gap equation reads 

a«/j(p/)=- E E 

i—{x,y,z} 7<5 

x (o^tiPf - p' f )4s n (p'f)MPf)) p , (3) 

where f 7 s{Pf) = TJ2T n J d£ P Eys(Pf,e n ;€p), F is the 
anomalous propagator, e n fermionic Matsubara fre- 
quencies, and e c is the usual frequency cut-off. The 
anisotropic interaction in Eq. {^) breaks spin rotational 
symmetry, but since each py-state is doubly Kramers- 
degenerate in zero-held, we can still decompose A Q ^ 
and F 7< 5 into (pseudo-)spin singlet (s) and (pseudo-)spin 
triplet (t) components [fl9f and arrive, for x = s, t, at 

A,( P/ ) = -{V x (p f - P'f) n{p' f ) h(p'f)) n , ( 4 ) 

with v. = x x + r/2, v t a = -x 1 - + r/2, v/- = -r/2. 

Isotropic spin fluctuations (x x = X z ) support triplet su- 
perconductivity only for nearly ferromagnetic enhance- 
ment. In the case of extreme anisotropy the coupling 
functions for singlet and triplet pairing with d-vector in 
z-dircction have equal sign and magnitude, V s — V t z — 
X z /2. In addition, a second triplet-pairing state with 
d _L z is possible. It couples via V f x = — x z /2. 

In order to study the symmetry of the superconducting 
state near T c q we determine numerically the eigenvalue 
spectrum of the integral kernel in Eq. ^, following Ref. 

. The resulting complete orthogonal set of basis func- 
tions yjj,(Pf) can be classified according to the irreducible 
representations (r) of the crystal group D^h- The corre- 
sponding eigenvalues, A^, determine the coupling con- 
stants for the /ith basis function in the symmetry chan- 
nel (r) plj . The most attractive (negative) eigenvalue in 
representation (r), A r = min /1 (A^) is eliminated in favor 
of a transition temperature for order parameter symme- 
try (r) in the usual way, T c y — 1.13e c exp(— l/|A r |). 

The dominant coupling constant A r ° determines the 
superconducting transition temperature T c q and the sym- 
metry (rg) of the superconducting phase near T c q. Once 




FIG. 1. The most attractive eigenvalues, A r , for the differ- 
ent irreducible D^h crystal symmetry group representations, 
(r), as a function of the magnetic incommensuration param- 
eter S at £ s fl = 4a. Even-parity representations are Ai g , 
A2 g , Big, and £?2 9 - The odd-parity states have either d \\ z 
(E2u(z)) or d _L z (E2u{xy)). The regions, separated by ver- 
tical lines, are marked by the representation of highest T c . 



an order parameter component in a certain representa- 
tion (r) nucleates, the other components fi(T) in the 
same representation are usually induced by the presence 
of the first component. Thus, physical transitions be- 
tween different superconducting phases only occur when 
additional symmetries are spontaneously broken. Such 
subdominant transitions are suppressed below the value 
T c r in the presence of a dominant order parameter. In 
the following we study which of the different possible 
phases nucleates first at given £ s a and 5 and determine 
points of accidental degeneracy between two different or- 
der parameter phases as a function of these parameters. 
In Fig. [j] we show the dependence of the attractive eigen- 
values on S for x = 0, fixing £ s a at 4.0a, a value which 
should closely correspond to the actual value in Sr2Ru04 
P|. The relevant phases near T c q are the even-parity 
one-dimensional irreducible representations Ai g and B\ g , 
which give spin-singlet superconductivity, and the odd- 
parity two-dimensional irreducible representation, E% u , 
rendering a spin-triplet channel for superconductivity. At 
small values of incommensuration, up to S ~ 0.35, the 
system prefers the -Bi ff -pairing channel. Above this point 
there is a region, 0.35 < S < 0.42, of spin-triplet pairing 
in the Ei u channel with d || z. Beyond 5 « 0.42, the pair- 
ing state is spin-singlet with A\ g symmetry, bounded by 
a narrow region of Bi ff -pairing starting around 8 ~ 0.61. 
In the large-(5 range there is again triplet pairing, but 
now with d _L z. The two accidental degeneracy points 
of interest for us are B\ g © E 2u for an incommensuration 
5 « 0.35 and A lg © E 2u for 6 w 0.42. Both points are 
remarkably close to the experimental value of 5 w 0.4 
and the theoretically predicted value of 6 w 0.33. Cal- 
culations with an additional component x^ resulted into 
accidental degeneracies even closer to each other. We 
also performed calculations for 5 w 0.4 and varying £ s a 
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showing that the presence of the accidental degeneracies 
near 8 = 0.4 is a robust feature for £ s g > 2a. 

In the following we study the possible low temperature 
superconducting phases close to the accidental degenera- 
cies. We concentrate on the three cases where either of 
the three symmetry channels is slightly dominant. We 
chose 8 — 0.3 for a dominant B\ g channel, 8 = 0.4 for a 
dominant E^ u channel, and 8 — 0.45 for a dominant A\ g 
channel. To determine the superconducting state at low 
temperatures we solve the non-linear gap equation (^). 
Expanding the order parameter with respect to the set 
of basis functions we obtain the order parameter 

components Aj^ for each representation. In the case of 
a mixed parity state a mixture of even parity, A' e )(py), 
and odd-parity, A^(pj), basis functions occurs as 

a ± ( P/ ) = y, ^y^ipi) ± E A (5) 

Once parity is broken, each of the doublets, ip\ip± and 
ifoipf, acquire separate order parameters, ip^ipi with 
A+(py) and r^ifo with A_(pA and A(py) has the re- 
quired anti-symmetry since A±(py) = — A=p(— py). As 
there may be several possible superconducting states, 
each state being a local minima in the free energy, we 
compute the free energy of each candidate state using 
the Serene- Rainer free energy functional p2[ , and select 
the state of lowest energy as the low-tcmpcrature phase. 
As the temperature evolution of the order parameter de- 
pends on T c ,r/T c0 = exp(l/A r — 1/A r °), it is dependent 
on the value of the most negative eigenvalue A r ° , which 
is determined by the parameter g 2 XQ- We have chosen 
<7 2 Xq/(27t) 2 = 1, leading to eigenvalues A r < 1, appro- 
priate for weak coupling. 

In Fig. H we plot the magnitude of the superconduct- 
ing gap as a function of position on the Fermi surface, 
A + (py)|, that minimizes the free energy at 8 = 0.4 and 
with £ s fl = 4a. We find three local minima: one ground 
state (G) and two metastable states (Ml and M2). The 
free energy difference at zero temperature between Ml 
and G corresponds to only 9% of the ground state con- 
densation energy (it amounts to 22% for M2). G is of 
pure £^2u-symmetry. Ml is symmetric around the p x and 
p y axis and has points with small gap values on the a- 
sheet. M2 is symmetric around the diagonals p x ± Py 
and has nodes on the a sheet. Also shown in Fig. 
is the total (angle-averaged) density of states (DOS) at 
the Fermi surface. The DOS is fully gapped for the 
ground state. The first metastable state shows a DOS 
with a much smaller excitation gap, and the DOS for the 
second metastable state shows low energy nodal excita- 
tions originating from the a sheet. The metastable states 
break spin-rotation symmetry and parity. The ampli- 
tudes A/f' ' have a zero relative phase for basis functions 
within the same parity, and a relative phase difference of 
±7r/2 for basis functions with opposite parity; thus, 




at zero temperature for the superconducting ground state 
(G) and the two metastable states (Ml, M2), calculated for 
8 = 0.4 and £ s fl = 4a. The three Fermi sheets (a, (3 and 7) 
are marked by thin white lines and |A+(jjy)| by the thickness 
of the black lines. The ground state is of pure E2u symme- 
try, the metastable states mix E% u , Ai g , and Bi g , resulting 
in mixed parity states. The corresponding total density of 
states, N(e), for each state is shown in the lower right panel. 

A_(py) = A + (pj)*. In the ground state the odd-parity 
basis functions with same eigenvalue have a 7r/2 relative 
phase difference, giving a l p x ±ip y '-state. We emphasize: 
all three solutions break time-reversal symmetry. 

Going to a 8 where superconductivity nucleates first in 
either the B\ g or the A\ g channel, we find that already 
the ground states are mixed parity states. The lowest en- 
ergy state is a nodeless time reversal symmetry breaking 
state with relative phases similar to the Ml state dis- 
cussed above. For both 8 = 0.3 and 8 — 0.45 we find 
metastable states, with nodes on all three Fermi surface 
sheets, that are of slightly higher free energies (9.0% and 
15.6% respectively) than the nodeless ground states. In 
Fig. H we show |A+(p*)| and the DOS for both states at 
8 = 0.3 and 8 = 0.45. For the metastable states the odd- 
parity components can develop higher order nodes, and 
remarkably, even whole arcs with vanishing or very small 
order parameter magnitude at the positions of the nodes 
in the even-parity components. At 8 = 0.3, for instance, 
the odd parity component has nodes along a diagonal in 
the Brillouin zone, coinciding with the d-wave nodes of 
the even parity component. In this case such arcs of al- 
most zero A(°\pf) occur on the 7-sheet, extending all 
over the two quadrants which contain the nodes of the 
total order parameter. 

The two components of the odd-parity order parameter 
transform like a vector in momentum space. Its direction 
is determined by the anisotropy introduced by a) the nor- 
mal state DOS and b) the presence of the even-parity 
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FIG. 3. The same as in Fig. g| for S = 0.3 (left) and 
for 8 — 0.45 (right). Top to bottom: gap magnitude of the 
metastable state (M), gap magnitude of the ground state (G), 
and the density of states for each corresponding state. 

Big or Ai g components. The mixed parity solutions for 
each incommensurability correspond to the alignment of 
this vector with high symmetry directions in the Bril- 
louin zone. Because the energy difference between the 
ground state and the state with nodes is only a fraction 
of the ground state condensation energy, fluctuations in 
the direction of this vector can be responsible for the 
presence of nodal excitations. Analyzing data of specific 
heat measurements || and of thermal conductivity ||, 
prompts that the Si'2Ru04 order parameter should have 
nodes This conclusion is further strengthened by re- 
cent measurements of the penetration depth showing a 
non-exponential low-temperature behavior . Based on 
the presented calculations, this implies that a mixed par- 
ity superconducting state, breaking spin-rotation symme- 
try, may be a candidate state for Sr2Ru04, and maybe 
even stabilized by additional interactions in the Hamil- 
tonian not considered here. 

In conclusion, we have shown that in the parameter 
region supported by experimental results even parity sin- 
glet and odd parity triplet pairing compete to have the 
highest transition temperature. Given the spin-spin cor- 
relation length of several lattice constants, as suggested 
by experiment, for an incommensuration near S = 0.35 
superconductivity nucleates in an accidentally degener- 



A\ g © E 2u - Both states are of mixed parity and break 
the T)±h crystal symmetry and time reversal symmetry. 
Within our model we find as ground states time rever- 
sal symmetry breaking nodeless states. However, close 
in free energy there exist metastable states of the order 
parameter with nodes, which may be stabilized by addi- 
tional interactions. 
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